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ABSTRACT 

Among its great findings, the IRAS mission showed the existence of an unidentified 
mid-IR feature around 21 /um. Since its discovery, this feature has been detected in ah 
C-rich proto-PNe of intermediate spectral type (A-G) and - weakly - in a few PNe and 
AGB stars, but the nature of its carriers remains unknown. In this paper, we show the 
detection of this feature in the spectra of three new stars transiting from the AGB to 
the FN stage obtained with the Spitzer Space Telescope. Following a recent suggestion, 
we try to model the SEDs of our targets with amorphous carbon and FeO, which might 
be responsible for the unidentified feature. The fit thus obtained is not completely 
satisfactory, since the shape of the feature is not well matched. In the attempt to relate 
the unidentified feature to other dust features, we retrieved mid-IR spectra of all the 
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21-fj,m sources currently known from ISO and Spitzer on-line archives and noticed a 
correlation between the flux emitted in the 21-//m feature and that emitted at 7 and 11 
Hm (PAH bands and HAC broad emission) . Such a correlation may point to a common 
nature of the carriers. 

Subject headings: Stars: AGB and post-AGB, circumstellar matter; Infrared: stars 



1. Introduction 

One of the major results of the Infrared Astronomical Satellite (IRAS) mission came from 
its Low Resolution Spectrometer (LRS) with the detection of a new broad feature around 21 



(jKwok et al. Ill989l ) in the spectra of stars that have gone through the Asymptotic Giant Branch 
(AGB) phase. The existence of this feature was not recognized at once, because it is weak in the 
first object in which it was observed, IRAS 22272+5435 , which was listed in the LRS Atlas as having 



a low-temperature continuum and silicate absorption (jVolk et al. 111993 ). In the first paper about 
this feature, four sources are listed (namely IRAS 07134+1005, 23304+6147, 22272+5435, and 
04296+3429); subsequent ground-based (United Kingdom InfraRed Telescope), airborne (Kuiper 
Airborne Observatory), and space-based (Infrared Space Observatory) observations confirmed those 
detections and also added new members to the group of 21-/im emitters. Observations with the 
Infrared Space Ob servatory (ISO) led to conclude that the feature actually peaks around 20.1 ;um 



(jVolk et al. Ill999l ). although it is still today commonly referred to as the "21 ^m feature". 



The feature has been detected almost exclusively in C-rich proto-Planetary Nebulae (PPNe), 
intermediate-mass stars (1-8 Mq) transiting from the AGB to the Planetary Nebula (PN) stage. 

During the AGB phase, stars lose most of their initial mass, with mass loss rates ranging from 
10^'' up to 10~^ Mq yr~^. The AGB stellar atmospheres are subject to periodic He-flashes, which 
imply that enriched material produced by nucleosynthesis in the bottom layers is carried to the 
surface. If this dredge-up of the star is efficient enough, the chemistry of the cooled material in the 
circumstellar envelope (CSE) will be based on C rather than O, as normally expected. Depending 
on their [C]/[0] ratio, we can then distinguish between 0-rich and C-rich CSEs. 

When the star has evolved past the tip of the AGB, the intense mass loss that characterizes 
the AGB has ceased, but the dust and gas in the GSE is typically so thick that the star cannot 
be detected at optical wavelengths. The expansion and subsequent dilution of the envelope leads 
to the optical detection of a post-AGB star, which typically has a double-peaked Spectral Energy 
Distribution (SED), as the result of the radiation from the central source and its CSE. 

The central star then evolves towards hotter temperatures, while its CSE continues to expand 
and dilute. If the central star does not evolve too slowly, a few thousand years after the AGB 
phase it is hot enough to ionize its CSE (T^ ^20-30x10^ K), which can lead to the destruction 
of molecules and even dust grains. The point at which the ionization of the CSE occurs is in fact 
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considered as the beginning of the Planetary Nebula phase. 



All of the objects showing the 21 nm feature a re C-rich and metal poor with enhanced abun- 
dances of s-processed elements (jHrivnak et al. and typically with F-G spectral type. Besides 
the unidentified feature, their mid-IR spectra typi cally show a broad feature around 30 /um that 



is usually attributed to MgS (jMolster et al. I l201Cl ) and has been found to consist of two features 



cente red at 26 and 33 ^m, though such a distinction is not observed in all sources (jHrivnak et al. 
2OO9I ). 



As of this writing, 16 PPNe are known to exhibit the unidentified 21-^m feature (iHrivnak et al. 



2001 



2009 ) . Besides these s ources, the feature h as been detected weakly in three PNe (jHony et al 



Volk et al. II2OO3I) and three A GB stars (|Volk et al. Il200d ). It has not been observed in IRAS 



01005+7910 (jCerrigone et al. II2009I ). a C-rich PPN with a ho t central star (BO I) that has already 
started to ionize its circumstellar shell (jCerrigone et al. II2OIOI ). Other hot post- AGB stars showing 
features of Polycyclic Aromatic Hydrocarbon s (PAH) in their mid-I R spectra and therefore likely 
to be C-rich do not exhibit the feature either (jCerrigone et al. II2009I ). This leads to the conclusion 
that the carrier must form during the AGB phase or very early afterwards, and is probably easily 
destroyed as the radiation field of the central star hardens. 

Since its first detection, many different molecules have been considered as possible carriers of 
the 21-//m feature. The connection with the C-rich nature of the envelopes makes it obvious to take 
into account C-bearing molecules, for instance TiC, doped-SiC, PAH, Hydrogenated Amorphous 
Carbon (HAG), and fullerenes. Yet, non C-bearing mole cules have been ind icated as possible 
carriers as well, for example SiS2, FeO, Fe203, and Fe304 (jZhang et al. II2009I ). The main issues 
that arise with almost all of the proposed carriers are the presence of secondary features that are 
not detected in the astrophysical spectra, or the anomalous abundances necessary to account for 
the observed intensity of the feature. 

In this paper, we present observations obtained within a program aimed at characterizing the 
dust emission in objects transiting fr om the AGB to t h e PN stage with A-G spectral types, as a 
follow-up to our project presented in ICerrigone et al. I (j2009l ). which focused on sources with hot 
central stars (B spectral type). While the whole sample of transition targets will be presented in 
a separate paper along with data from other Spitzer programs (Hart et al, in preparation), here 
we will focus on three new sources that show the unidentified 21-/im feature. These are the only 
targets in our sample, where we find evidence for this feature. 



Observations 



Th e new observations presented in this paper were carried o ut with the Infrared Array Camera 
(IRAC: lFazio et al. II2OO4I) and the InfraRed Spectrograph (IRS: lHouck et al. II2OO4I ) on-board the 
Spitzer Space Telescope (jWerner et al. II2004I ). within Program 50116 (PI: G. Fazio) in March and 
April 2009. The IRS had four modules with different spectral resolutions. The Short-Low (SL) and 
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Long-Low (LL) modules covered the 5.2-14.5 and 14.0-38.0 fim ranges respectively, with spectral 
resolution R~60-120. The Short-High (SH) and Long-High (LH) modules covered the 9.9-19.6 
and 18.7-37.2 //m ranges with R~600. Our spectral data were typically acquired in low resolution 
mode. Targets with IRAS fluxes that would imply saturation with the IRS at longer wavelengths 
were observed in low resolution only between 5 and 14 (SL modules), while the 10-38 //m range 
was observed in high spectral resolution (SH and LH modules). In the latter case, a dedicated 
background observation was also performed for each target. For the three sources presented here, 
low resolution observations were performed for IRAS 13245-5036 and 15482-5741, while low and 
high resolution data were collected for IRAS 14429-4539. IRAS 13245-5036 was not centered on 
the SL slit, therefore only LL data are available for this target. 

Basic Calibrated Data were retrieved from the Spitzer web archive. The observations of IRAS 
13245-5036 were processed with the 18.7.0 pipeline by the Spitzer Science Center (SSC), while 
version 18.18.0 was applied to 14429-4539 and 15482-5741. The data were reduced following the 
standard recipe in the on-line IRS Cookbook: the background subtraction was performed from the 
off-target observations for low resolution data (i.e., when one module is on the target, the other one 
is observing off-source) and from the dedicated background observation for the high resolution data 
sets. The spectra w ere cleaned using IRSCLEAN, extracted and averaged together in SMARlQ 



(jHigdon et al. II2004I ). The hi gh spectral resolution dat a were extracted over the full slit, while the 



optimal extraction algorithm ([Lebouteiller et al. Il2010l ) was applied to the SL and LL observations. 



The IRAC observations were performed in high dynamic range mode with a small dither 
pattern, as very extended emission was not expected. In this observing mode, two sets of BCDs 
are created by default, corresponding to a long and a short exposure; in our case the exposures 
were respectively 12s and 0.6s long. All of the IRAC BCDs retrieved from the Spitzer archive were 
processed with the S18.18.0 pipeline. The BCDs were treated for artifact correction with the IDL 
code available from the SSC web site and then multi-frame aperture photometry was performed 
with the MOsaicker and Point source Extractor (MOPEX) using a 10-pixel aperture and 12-20 pixel 
annulus, which is the standard for the IRAC calibration. As the IRAC absolute flux calibrators 
are known to about 2-3%, and IRAC repeatability is at about the 1.5% level, we assume in the 
following 5% as our photometric error. Since our targets are quite bright for IRAC, the photometry 
extraction was performed on the short frames, to avoid any saturation effect, which was evident in 
IRAS 14429-4539 at 8.0 /im even in the short frames. The flux densities are listed in Table [H 



^ SMART was developed by the IRS Team at Cornell University and is available through the Spitzer Science Center 
at Caltech. 
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3. Spectra and SED modeling 



Figure [T] shows the results of our spectral observations. For IRAS 13245-5036, only long- 
wavelength data are shown, because the short- wavelength slit (perpendicular to the long-wavelength 
one) was off the target. The 21-/im feature is evident in all of the three targets. IRAS 14429-4539 
and 15482-5741 also clearly show the mid-IR features attributed to PAHs, although the balance 
among the features in the PAH bands seems to be very different in the two sources. In IRAS 14429- 
4539 these features resemble closely those detected in PNe and other C-rich stars, while usually 
this is not the case for the 21-/im sources, where the PAH features in the 10-17-/im range (due to 
CH out-of-plane bending modes) are often blended together in a plateau (see for example IRAS 
22223-1-4327 in Figure [6|) and are as strong or even stroi iger than t he fea tures in the 5-10-/im range 
(due to CC stretch and CH in-plane bending modes; iTielens. I (j2008l )). as for example in IRAS 
22574-1-6609. Figure [2] shows the normalized profiles of the 21-//m features in our targets and in 
IRAS 23304-1-6147, which can be considered as prototypical. As usually observed, the feature has 
a relatively regular spectral shape, with no hints for substructure. 

To estimate the relative contributions of the PAH features within the emission bands, we have 



fitted the spectra using the publicly- available IDL tool PAHfit (jSmith Sz Drain 1120071 ). specifically 
created to fit PAH features in low-resolution Spitzer spectra. The spectra were binned to a step 
size of 0.169 /im. 

PAHfit performs a fit to the spectra assuming the continuum is given by a combination of gray 
bodies with set temperatures of 300, 200, 135, 90, 65, 50, 40, and 35 K and a stellar continuum 
given by a Planck curve at a temperature input by the user. The method of fitting gray bodies is 
based on the fact that at these wavelengths the underlying continuum is a combination of radiation 
from the central star and the dust surrounding it. The dust layer extends outwards as far as it 
cannot be distinguished from the ISM, therefore it is characterized by a continuum of temperatures, 
typically in the range of a few 100 K to 10-20 K. In our two cases, the software used gray bodies 
at 300, 200, and 135 K for IRAS 14429-4539 and at 300, 135, 90, and 35 K for 15482-5741. The 
fact that PAHfit returns different temperature combinations for the two sources should not be 
over-interpreted, as we are still fitting only a si nall spectral r a nge. For the central stars, since all 
of our targets have been optically classified in ISuarez et al. I (j2006l ). we adopt the temperatures 
approximately corresponding to their spectral classifications: 7000 K for 14429-4539 and 6500 K 



Target 


3.6 fim 


4.5 fiiR 


5.8 ^m 


8.0 fim 




mJy 


mJy 


mJy 


mJy 


13245-5036 


21.08 


17.24 


71.21 


440.9 


14429-4539 


573.3 


1132 


3256 




15482-5741 


40.41 


32.37 


47.39 


209.3 



Table 1: Photometry of our targets obtained with IRAC. The absolute accuracy is about 5%. 
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for 15482-5741. We list the output obtained from PAHfit in Tables [2] and [3] and the fits are shown 
in Figure [3l Besides PAH features, PAHfit also tries to fit several atomic lines, which we have 
included in our tables. These lines are blended with the PAH features: their presence is supposed 
to obtain a good fit. In IRAS 15482-5741, two atomic lines are particularly weak and marginal: no 
error estimates are given for these lines. 

As a further tool to analyze the properties o f the envelopes in these stars, we performed 



radiative transfer modeling using the code DUSTY (jlvezic et al. Ill999l ). which solves the equation 



of radiation transfer in a spherical shell of dust grains surrounding a central star. For the grain 
sizes, we adopted a standard MRN (jMathis et al. 1119771 ) distribution of radii varying as a~^'^ with 
amin = 0.005 /im and amax = 0.25 /im, a being the radius of each grain (assumed spherical). A 
radial density distribution going as was implemented in all models and Planck curves were used 
to model the radiation from the central stars. The choice of a density distribution varying as r~^ 
is based on the expectation of a constant mass-loss rate. From the equation of mass continuity, we 
have that 

M 



p{r) 



(1) 



47rr2y^ 

where Voq is the wind terminal velocity. As a consequence, if the mass loss process is steady, the 
CSE will be characterized by p{r) ~ r^^. Towards the end of the AGB, the star is supposed to 
undergo a short (a few 10'^ yr) period of enhanced mass loss called superwind phase, therefore 
departures from an r~^ distribution may be found in these s tars. Once a sat i sfacto ry fit to the data 
was achieved, we reddened the modelled SED according to Cardelli et al. (jl989l ). taking Ay as a 
free parameter. We also summed to the DUSTY model, an empirical fit to the 30-^m MgS feature 
obtained with an asymmetrical Gaussian. This was done to better reproduce the e mission under 



the 21 -j^m feature, which we tried to model as due to FeO, as recently suggested bv IZhang et al. 



(j2009l ). DUSTY handles mixtures of chemical species as simulating a single- type grain constructed 
from an average. For both amorphous carbon and FeO, we used the default sets of optical constants 
that come with DUSTY. The details of the modeling are summarized in Table H] and the models 
are shown in Figure HI 



The observation al points us ed in modeling are fro m our IR AC photometry and th e GSC2.2 



Lasker et al. 



200^1, NOMAD dZacharias et al. Ibood), D ENIS dEpchtein et al. Ill997ll. 2MASS 



Skrutskie et al. 



20061 ). AKARI ( 



Murakami et al. 



20071 1. and IRAS (jNeugebauer et al. I \l9m 



catalogs. 

As often observed in post- AGB stars, we found that while the far-IR observational points can 
be fitted well, the model typically underestimates the radiation at near- and mid-IR wavelengths. 
This implies that a larger amount of hot dust is present in the CSE than what can be found in a 
smooth spherical shell. To account for such an emission component, we have tried different density 
distributions (r^^ and r^^), but found that a satisfactory fit can be obtained by running the code 



twice with p{r) 



-2 



In the first run we try to reproduce the data at wavelengths below about 



10 /im, then in the second run we use as a central source the SED obtained as the output of 
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Feature Wavelength Peak intensity Flux 

urn Jy 10-1^ W m-2 



PAH features 






5.3 


5.27 


0.24 ± 0.02 


7.4 ± 0.5 


6.2 


6.22 


3.84 ± 0.06 


87 ±1 


6.7 


6.69 


0.515 ± 0.009 


25.4 ± 0.4 


7.4 


7.42 


0.71 ± 0.03 


56 ±2 


7.8 


7.85 


8.16 ±0.06 


260 ±2 


8.3 


8.33 


1.40 ±0.02 


39.4 ±0.5 


8.6 


8.61 


4.32 ± 0.02 


92.2 ±0.4 


10.7 


10.68 


1.48 ±0.02 


13.0 ±0.2 


11.3 


11.33 


8.31 ±0.07 


111 ± 1 


12.0 


11.99 


3.75 ± 0.04 


66.3 ±0.6 


12.6 


12.62 


3.75 ± 0.04 


58.7 ±0.6 


13.5 


13.48 


1.67 ±0.04 


23.3 ±0.6 


14.2 


14.19 


0.86 ± 0.07 


7.1 ±0.6 


17.4 


17.37 


0.3 ±0.1 


0.9 ±0.4 


Atomic lines 








H2 S(5) 


6.90 


1.09 ±0.07 


4.24 ±0.27 


[Ar II] 


7.01 


0.73 ± 0.04 


2.77 ±0.17 


H2 S(4) 


8.10 


3.66 ± 0.09 


19.8 ±0.5 


[Ar III] 


8.94 


0.93 ± 0.02 


4.06 ± 0.09 


[S IV] 


10.53 


0.54 ±0.03 


1.71 ±0.10 


H2 S(2) 


12.30 


0.94 ±0.05 


2.18 ±0.12 


[Ne II] 


12.86 


0.41 ± 0.06 


0.88 ±0.13 


H2 S(l) 


17.08 


0.90 ±0.14 


1.52 ±0.24 



Table 2: Output obtained from the IDL tool PAHfit for IRAS 14429-4539 
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Feature 


Wavelength 


Peak intensity 


r lux 




fim 


Tit 


iu w m 


PAH features 






5.3 


5.27 


0.008 ± 0.002 


0.25 ± 0.07 


6.2 


6.22 


0.053 ± 0.002 


1.21 ±0.04 


6.7 


6.69 


0.029 ± 0.001 


1.42 ± 0.06 


7.4 


7.42 


0.090 ± 0.002 


7.2 ±0.1 


7.8 


7.85 


0.110 ±0.002 


3.51 ±0.07 


8.3 


8.33 


0.106 ±0.002 


3.00 ± 0.05 


8.6 


8.61 


0.066 ± 0.001 


1.40 ± 0.03 


10.7 


10.68 


0.009 ± 0.004 


0.08 ± 0.04 


11.3 


11.33 


0.185 ±0.005 


2.46 ± 0.06 


12.0 


11.99 


0.128 ± 0.004 


2.27 ±0.07 


12.6 


12.62 


0.110 ±0.006 


1.72 ±0.09 


13.5 


13.48 


0.055 ± 0.004 


0.77 ±0.05 


14.2 


14.19 


0.071 ± 0.007 


0.59 ± 0.06 


17.4 


17.37 


0.01 ±0.01 


0.045 ± 0.04 


Atomic lines 






H2 S(5) 


6.916 


0.0648 ± n.a 


0.22 ± n.a. 


[Ar II] 


7.002 


0.0336 ± n.a. 


0.11 ± n.a. 


H2 S(4) 


8.076 


0.073 ± 0.003 


0.39 ± 0.02 


[Ar III] 


8.994 


0.060 ± 0.003 


0.26 ± 0.01 


H2 S(2) 


12.264 


0.06 ± 0.01 


0.15 ±0.02 


[Nc III] 


15.505 


0.06 ± 0.02 


0.12 ±0.04 


H2 S(l) 


17.085 


0.011 ±0.003 


0.018 ± 0.005 



Table 3: Output obtained from the IDL tool PAHfit for IRAS 15482-5741. 



Target 


T^ 


[FeO] 


L 
d 


rnin 
-^dust 


Rout 
Rin 


TV 


Av 




K 


% 


L0 kpc^^ 


K 








13245-5036 


8500 


9 




220 


1.1 


0.25 








9 


267 


100 


1.1 


0.25 


1.05 


14429-4539 


7000 


7 




500 


10 


0.6 








7 


673 


160 


20 


1.6 


4.1 


15482-5741 


6500 


20 




255 


10 


0.07 








20 


267 


120 


30 


0.55 


2.3 



Tabic 4: Details of the modeling performed with DUSTY. Except abundances, the same parameter 
values were used to generate models with and without FeO. 
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the first run, which contains the central star and the hot dust. Such a density picture resembles 
a circumstellar environment where the mass-loss process is not constant with time, but subject 
to episodic enhancements. The presence of a hot- dust component can als o be interpreted as an 
indication for a circumstellar (circumbinary?) disk (jPe Ruvter et al. 1120061 ). 



Our attempt to reproduce the 21-/im feature with FeO following IZhang et al. I (|2009l ) is not 



satisfactory. In particular, the spectral shape is not well reproduced and large abundances (~10% 
in mass) of FeO are necessary. Software allowing for a continuous range of ellipsoi dal grain shapes 



has b een found to account better for spectral shapes of dust features in general ([Kemper et al 
2OO2I ). It is possible that a similar approach may improve the match between the model and the 



observational data also in our targets. Currently, this is not possible with DUSTY. 



4. A link between the 21-//ni feature and hydrocarbon emission? 



It is known that the 21-/xm feature is often - but not always - observed along with mid-IR PAH 



bands (IVolk et al. Ill999l ). Another source of mid-IR emission in the form of a very broad feature 



resembling an underlying continuum is given by HAG, whose molecular chains are considered as 
the precursors of PAHs. Given this qualitative observational link between these features, we have 
searched for a quantitative correlation. 

We have retrieved from the ISO and Spitzer archives all of the available IR spectra approxi- 
mately covering the 5-40 nm range for the so-far known 21-/xm PPNe. Table [5] lists the data found 
in the archives and used in the following analysis. Data for IRAS 05113-1-1347, 06530-0213, and 
07430-1-1115 were not used because they do not cover the whol e range from 5 /im, w hile the data set 
for IRAS 19477-1-2401 was not used because of its low quality ([Hrivnak et al. II2OOOI ). When Spitzer 
data were available, we carried out our analysis on these, after combining them with the ISO data 
below 5 Atm. Post-BGD data were retrieved from the Spitzer archive, while the highly-processed 
data bv lFrieswiik et al. I (|2007l ) were retrieved from the ISO web archive. Glear outliers in the ISO 
data were manually flagged out in SMART, then a standard cj-clip run was performed and the 
whole procedure (manual flagging and cj-clipping) was repeated a second time, before averaging 
the data together. For all of the spectra the same spectral step of 0.169 /im was set in SMART. 
The spectra thus obtained are shown in Figures [5] and El 

To enhance the dust features, we fitted the underlying continuum as a linear combination of 
four gray bodies w ith emissivity index /? = 2. In doing so, we used the IDL routine MPFITFUN 
(jMarkwardt. II2008I ). which fits a user-defined function to the data. Although four gray bodies were 
always included in the fitting procedure, in most cases the best-fitting combination had a gray 
body multiplied by such a low factor that made it negligible, thus obtaining a best-fit curve with 
three gray bodies only. The same behavior was observed in trials with more than four gray-body 
components: only three or four of them were significant. As mentioned previously, the infrared 
continuum in post-AGB stars is due to a combination of stellar radiation and thermal emission from 
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Target 


Data sets 


Range (/im) 


Z02229+6208 


SWS^ 


l- 


-45 


04296+3429 


SL^ , SH^ , LH^ 


5~ 


-37 


05113+1347 


SH, LH 


10 


-37 


05341+0852 


SWS, SH, LH 


2- 


-45 


06530-0213 


SH, LH 


10 


-37 


07134+1005 


SWS 


2- 


-45 


07430+1115 


SH, LH 


10 


-37 


16594-4656 


SWS 


z- 


-40 


19477+2401 


SWS 


2- 


-45 


19500-1709 


SWS 


2- 


-45 


20000+3239 


SWS 


2- 


-45 


AFGL 2688 


SWS 


2- 


-45 


22223+4327 


SL, SH, LH 


5- 


-37 


22272+5435 


SWS 


2- 


-45 


22574+6609 


SWS, SH, LH 


2- 


-45 


23304+6147 


SWS, SH, LH 


2- 


-45 



Table 5: Summary of the archive observations available for all of the known 21-;um sources besides 
the new detections reported in this work. Targets not included in our analysis are in italics. 



°- Short Wavelength Spectrometer (ISO) 
^ Short Low IRS 
" Short High IRS 
Long High IRS 

dust, with the latter typically much stronger than the former. To avoid physically unreasonable 
solutions in the fit procedure, we have therefore limited the range of values investigated by the 
fitting routine to < T < 15000 K, which takes into account temperature values appropriate for 
both the dust (dust grains sublimate around 1000-2000 K) and the central star. The temperatures 
thus obtained are listed in Table [6] and the fitting curves are shown in Figures [5] and [6j We find that 
in all sources the emission is dominated by gray bodies reproducing dust emission (because of their 
temperatures) and stellar contributions are negligible. The continuum obtained as a combination of 
gray bodies is typically unable to properly isolate the feature at 21 //m, because of the substantial 
contributions from the features in the 26-33 /um range, which are often blended together. 

Emission in the 5.5-18 //m range is present in all of the targets, as shown in the left-side 
close-ups in Figures [5] and [6l After converting the spectra into W m"^ /im"^, we integrated the 
continuum-subtracted spectra to estimate the flux in the 5.5-18 fim and 18.5-23.5 //m ranges (the 
exact edges of the ranges were adapted to the emission in each source, turning out in differences of 
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about 0.3 fj,m). To assess the emission in the 21-nm feature, a straight basehne was fitted to its red 
and blue edges and then subtracted before performing the integration, as shown in two examples 
in Figure [71 The emission in the 5-18 //m range is mostly due to PAHs but contributions from 
HACs are possible, while the 18.5-23.5 //m range is obviously dominated by the 21-/nm feature. In 
Figure [51 we plot the flux estimated in the 5.5-18 //m bands, indicated as the sum of the PAH 7 
and 11 fim bands, versus that in the 21-/im feature. 

As can be seen in the plot, a correlation between the fluxes in the two wavelength ranges 
is evident; we can calculate a linear correlation coefficient of 0.79ib0.07. An explanation for the 
correlation may be that the abundances of the hydrocarbons and of the carrier of the 21-/im feature 
are linked together, possibly because of a common formation process and therefore the correlation 
would just reflect the different amounts of hydrocarbons and 21-/xm carriers produced in the CSEs. 

The colors of the dots in Figure [HI are proportional to the temperatures of the central stars, as 
derived from their spectral classifications, and range from red (coldest) to blue (hottest). No clear 
correlation is observed between the intensity of the features and the stellar temperature. 

In principle, if the observed correlation were due only to a common excitation process for the 
molecules involved and then essentially to the hardening of the radiation field, we would expect to 
find a correlation with the temperature of the central star. There may be several reasons for this 
lack of correlation. For example, extinction effects can play a role. Also, one of the explanations 
for the existence of the 21-//m feature implies radiation-induced decomposition of large HAC grains 



into smaller PAHs (jScott et al. II1997I ). In this scenario, the carrier of the 21-/xm feature would be a 
short-lived decomposition product. The efficiency of the processing of large grains would depend on 
both Tgff and the location of the dust layer. A massive star will be hot enough to cause processing 
earlier than a low-mass star, with its dust still close to the central star. A low- mass star will reach 
the same temperature at a later age, with its dust farther away and will therefore need an even 
larger temperature to attain comparable processing. 



5. Summary 

One of the main results of the IRAS mission was its discovery of the existence of a broad mid- 
IR feature around 20 /xm. The feature has been detected almost exclusively in the circumstellar 
envelopes of C-rich proto-PNe, with a few weak detections in their precursors (AGB stars) and 
successors (PNe). The nature of the feature is still unclear today, but it must be linked to the 
C-rich nature of the sources: its carriers are likely to be produced towards the very end of the AGB 
and are easily destroyed by UV radiation. 

In this paper we report about the detection of three new C-rich proto-PNe that exhibit this 
unidentified feature through observations performed with the IRS instrument on-board Spitzer. 
Besides the 21 /im feature, our targets also show mid-IR bands from PAHs and a broad feature 
around 30 fiui commonly attributed to MgS. As the dispersion of the dust formed during the AGB 
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creates a very extended shell around the central star, it can be expected that evidence for such 
large dust regions can be found in infrared images, but our IRAC observations do not point to any 
extended emission. 

By a combination of catalog data and our new observations, wc reconstruct the SEDs of our 
targets and model them with the code DUSTY. With a r^^ density distribution, a good match of 
the observational points is possible, if the model assumes a circumstellar environment where the 
mass loss process has not been constant with time. The use of FeO as carrier of the 21-/xm feature 
does not accurately reproduce the spectral shape of the feature. 

We have retrieved the infrared spectra of the other 21-//m sources available in the ISO and 
Spitzer archives. A gray-body continuum has been subtracted from the spectra to enhance the 
features. After continuum subtraction, we have calculated the flux in the features from 5 to 18 fim, 
due to hydrocarbons, and that in the 21-/xm feature, finding a clear correlation between them. Wc 
have not found any correlation with stellar temperatures. The correlation found between the 5-18 
and 21-/Lim fluxes can point to a common origin of the carriers of the features. 
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Facilities: ISO (SWS), Spitzer. 
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Target 






K 



Z02229+6208 


770 


280 


131 


68 


05341+0852 


639 


148 


57 




07134+1005 


1151 


198 


73 




16594-4656 


710 


159 


75 




19500-1709 


1495 


537 


125 


67 


20000+3239 


848 


321 


138 


64 


22272+5435 


970 


293 


142 


67 


22574+6609 


350 


130 


55 




23304+6147 


886 


208 


68 




AFGL 2688 


233 


117 


65 


40 


04296+3429 


387 


147 


60 




22223+4327 


251 


101 


55 




13245-5036 


509 


114 


67 




14429-4539 


355 


147 


60 




15482-5741 


498 


208 


154 


60 



Table 6: Temperatures of the gray bodies used to fit the underlying continuum. 
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Fig. 1. — The spectra obtained with Spitzer are shown on the right side, where a gray-body fit 
to the continuum (light-gray solid line) and the single gray bodies (light-gray dashed lines) are 
included. On the left side, a close-up of the 5-25 /xm range, after subtraction of the continuum. 
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Fig. 2. — Normalized profiles of the 21-/im feature in our targets compared with that in IRAS 
23304+6147. 



IRAS 14429-4539 IRAS 15482-5741 




Wavelength (jim) Wavelength (^tm) 

Fig. 3. — Fits to the PAH features obtained with the IDL tool PAHfit. Single gray bodies are 
shown in red and their total as a thick gray line. PAH features are shown in blue and atomic lines 
in purple. The sum of all components is shown as a green line. 
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Fig. 4. — DUSTY models obtained for our targets. The DUSTY output is shown as a dashed hue; 
a sohd hne representing the final model SED (DUSTY output plus an empirical fit to the 30 nm 
MgS feature) has been also plotted. On the left side, only amorphous carbon grains were used, 
whil e the models o r i the right side include a mixture of amorphous C and FeO. Optical spectra 
from lSuarez et al. I (|2006l ) and our Spitzer spectra are shown in green. 
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Fig. 5. — Infrared spectra of the 21-fj,m sources known. A gray-body continuum is shown as a gray 
sohd hne and the single gray bodies as dashed gray hues in the right-side plots. On the left side, a 
close up of the 3-24 fim range after continuum subtraction is shown. Gray-body temperatures are 
listed in Table O 
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Fig. 6. — As in Fig [5l 
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IRAS Z0222g+6208 IRAS 23304+6147 




Fig. 7. — Two examples of the straight basehne (dashed hne) subtracted from the 21-//m feature. 




Fig. 8. — Integrated flux emitted between 5 and 18 //m (containing the 7 and 11 fj,m PAH fea- 
tures), plotted versus the integrated flux in the 21-//m feature. The color coding depends on the 
temperature of the central star, with red being colder and blue hotter. 



